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摘 要 : 盐 渍 化 是 构成 绿洲 农业 低产 的 主要 原因 之 一 ,也 是 农业 
素 。 为 提高 盐 渍 化 耕地 生产 力 , 促 进 绿洲 农业 的 可 持续 发 展 ,以 喀什 噶 尔 绿洲 耕地 为 研究 对 象 , 利 


发 和 可 持续 发 展 的 重要 限制 因 


H| Landsat 8 OLI 遥 感 影像 数据 提取 还 感 指数 20 个 ,利用 土地 利用 数据 计算 研究 区 耕地 开垦 年 限 ， 
用 线性 拟 合 的 方法 将 用 植被 光合 作用 模型 (VPM) 模 拟 的 植被 净 初 级 生产 力 (NPP) 数 据 进 行 降 尺 
度 ,将 遥感 指数 同 土壤 采样 及 实测 数据 进行 相关 分 析 , 得 到 优选 的 遥感 特征 变量 ,再 用 粒子 群 优化 
算法 (PSO) 优 化 的 概率 神经 网 络 (PNN ) 模 型 进行 盐 渍 化 程度 分 类 ,得 到 研究 区 耕地 盐 渍 化 等 级 分 


布 情况 ,后 与 研究 区 耕地 


星 年 限 和 NPP 进 行 合 加 分 析 。 结 果 表 明 :(1) 选取 增强 型 植被 指数 


(EVI) 、 盐 分 指数 2(SL)、 湿 度 指数 (WI)、MSAVI-WI-SI 特 征 空间 (MWSI) 波段 6(B6,2.11~2.29 um) 
5 个 遥感 参量 通过 PSO-PNN 模型 进行 盐 渍 化 程度 反 演 准确 率 约 为 80%。(2) HMA BER MA BY 


区 域 盐 渍 化 程度 越 低 。 新 开垦 的 耕地 主要 分 布 在 研究 区 东部 ,而 研究 区 


部 大 都 为 开垦 年 限 在 45$ a 


以 上 的 老 绿洲 农业 区 。(3) 耕地 盐 渍 化 严重 降低 了 耕地 农作物 生产 力 。 研 究 区 耕地 NPP 较 高 的 区 
域 大 都 分 布 在 西部 , 较 低 的 区 域 大 都 分 布 在 东部 ,与 盐 渍 化 程度 等 级 分 布 大 致 相反 。 上 述 研 究 方 
法 与 结果 可 为 后 续 使 用 遥感 参量 进行 盐 渍 化 反 演 的 研究 提供 参考 ,对 干旱 半 干 旱 区 的 盐 渍 化 耕地 


改良 具有 一 定 的 参考 意义 。 
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文章 编号 : 


盐 渍 土 是 重要 的 耕地 后 备 资源 ,土壤 盐 渍 化 
和 次 生 盐 渍 化 使 人 们 赖 以 生存 的 耕地 遭 到 破 环 ,如 
果 不 及 时 阻止 耕地 盐 江 化 的 扩散 ,将 会 对 当地 经 济 
和 社会 发 展 造成 巨大 威胁 。 目 前 ,土壤 盐 清 化 问 
题 严重 影响 了 喀什 噶 尔 绿洲 地 区 农业 的 稳定 与 可 
持续 发 展 , 喀 什 噶 尔 绿洲 耕地 占 全 杜 总 灌溉 面积 的 
1/10, 亩 均 水 量 低 于 全 疆 平 均 水 平 一。 喀什 噶 尔 绿 
洲 由 于 气候 干旱 少雨 .地形 闭 塞 、 水 利 设施 不 配套 
等 原因 ,形成 了 盐 渍 化 较 重 . 盐 渍 化 耕地 面积 较 大 
的 绿洲 农业 区 ””。 对 大 面积 农田 盐分 信息 的 监测 
能 够 为 喀什 噶 尔 绿洲 开展 流域 规划 合理 配置 农业 
资源 提供 科学 依据 ,对 区 域 生 态 保护 和 农业 可 持续 
发 展 具有 重要 意义 。 


收 稿 日 期 : 2022-01-11; 修订 日 期 : 2022-02-11 


对 于 土壤 盐 汪 化 监测 ,传统 的 方法 主要 是 通过 
采集 土壤 样品 进行 室内 化 验 分 析 。 这 种 方法 流程 
复杂 、 耗 时 耗 力 ,无 法 进行 大 面积 动态 监测 。 而 目 
前 流行 的 遥感 技术 由 于 其 获取 简单 .限制 少 ,受到 
了 研究 人 员 的 青睐 , 越 来 越 多 的 人 利用 遥感 技术 
进行 盐 渍 化 信息 的 获取 7 ,并 不 断 取 得 了 一 些 进 
展 "”。 另 外 ,土壤 盐 渍 化 过 程 是 一 个 动态 的 ` 强 非 
线性 的 过 程 , 受 多 种 不 同 因素 “相互 作用 的 影响 ， 
这 给 如 今 的 盐 渍 化 评估 工作 带 来 了 困难 。 近 年 来 ， 
机 噩 学 习 方 法 由 于 其 能 够 自主 学 习 数据 处 理 速 度 
快 等 特点 被 越 来 越 多 地 应 用 于 盐 溃 化 研究 中 “”， 
提高 了 盐 溃 化 反 演 的 精度 ,同时 节约 了 大 量 人 力 物 
力 。 概 率 神经 网 络 (PNN) 是 一 种 新 型 非 线 性 神经 网 
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络 模 型 PNN 结构 简单 ,训练 时 间 短 ,需要 人 工 
调整 的 参数 少 ,判定 率 高 ,不 易 陷 人 局 部 极 值 ”， 
此 ,PNN 经 常 被 用 于 解决 分 类 问题 。 由 于 算法 简单 、 
搜索 速度 快 等 优点 ,粒子 群 优化 算法 (PSO) 被 广泛 
应 用 于 各 算法 优化 中 。 为 快速 摸 清 喀什 噶 尔 绿洲 耕 
地 盐 渍 化 状况 ,探究 盐 渍 化 与 植被 生产 力 和 耕地 开 
蜡 年 限 的 关系 ,本 文采 用 PSO 优化 的 PNN 模型 分 析 
喀什 噶 尔 绿洲 耕地 盐 渍 化 程度 分 布 , 以 期 对 当地 盐 
渍 化 耕地 管理 和 生态 系统 可 持续 发 展 提供 参考 。 


1 研究 区 概况 


研究 区 喀什 噶 尔 绿洲 (75"24'~78*00' 下 ,38"25'~ 
39?"56'N) 位 于 欧 亚 大 陆 内 部 ,新疆 西部 ,面积 4.62x 
10 kn? , 占 流 域 面积 4.02% , 占 新 疆 绿洲 面积 7.5%。 
喀什 噶 尔 绿洲 是 典型 的 暖 温带 大 陆 干 旱 气 候 , 降 水 
稀少 ,夏季 气候 干旱 炎热 .蒸发 强烈 ,冬季 寒冷 . 光 
照 充足 , 层 夜 温差 大 3。 喀 什 嘲 尔 绿洲 光照 资源 让 
富 , 适 宜 瓜 果 等 作物 的 生长 ,通过 农业 结构 的 调整 ， 
当地 农业 发 展 潜力 较 大 5 , 现 已 开发 利用 耕地 26.4x 
10'hm。 喀 什 噶 尔 绿洲 具有 悠久 的 农业 历史 , 主 产 
的 棉花 ,粮食 .特色 果品 在 全 疆 占 有 重要 地 位 。 


2 数据 与 方法 


2.1 数据 来 源 

2.1.1 ERYR 本 文采 用 的 遥感 影像 数据 源 
为 美国 陆地 资源 卫星 Landsat 系列 2020 4E. Landsat 8 
OLI 影 像 , 均 由 美国 地 质 调 查 局 USGS 网 站 (https:/ 
glovis.usgs.gov/) 下 载 获 得 ,获取 时 间 为 4 一 9 月 。 
2.1.2 耕地 开 圣 年 限 ”土地 利用 数据 来 自 中 国 科学 
院 资 源 环境 遥感 调查 建立 的 全 国土 地 利用 数据 集 
(NLUD) (https://www.resde.cn/) 。 利 用 1975—2020 
年 土地 利用 数据 空间 分 析 得 出 的 研究 区 耕地 开垦 
年 限 。 耕 地 开 星 年 限 具 体 计 算 方法 为 :通过 将 研究 
区 1975 一 2020 年 共 8 期 土地 利用 数据 的 徐 加 分 析 ， 
得 到 研究 区 耕地 开 星 年 限 ,并 分 为 5 个 时 段 ,分 别 为 
0-10 a, 10~20 a,20~30 a, 30~45 a %45a DJ Eo 

2.13 土壤 数据 采样 及 实验 测定 ”利用 GPS 定 位 在 
喀什 噶 尔 绿洲 耕地 共 选 取 130 个 代表 性 采样 点 (图 
1) ,并 在 2021 年 5 月 进行 土壤 采样 工作 ,130 个 采样 
点 的 布设 过 程 综合 考 虑 了 土壤 类 型 .土地 开 殿 年 限 
和 作物 类 型 ,其 中 ,棉花 地 82 个 .小麦 地 17 个 未 耕 


1135.4 
0 30 km 
l 


1 研究 区 采样 点 分 布 
Fig. 1 Distribution of sampling points in the study area 


种 地 7 个 .其 他 耕地 24 个 。 然 后 在 实验 室 对 130 个 
土 样 的 总 盐 RIRE OH) 八大 离子 等 指标 进行 了 
平行 化 验 分 析 ,测定 方法 均 采 用 国家 标准 方法 。 
2.1.4 NPP &4d& NPP 是 指 植被 在 单位 时 间 和 面积 
内 积累 的 有 机 二 物质 总 量 ,包括 植物 树枝 .叶子 和 
根部 以 及 植物 死亡 部 分 的 生产 中 。NPP 可 以 在 均匀 
尺度 上 反映 生态 系统 生产 力 ”“。 本 文采 用 中 国 科 
学 院 完 成 的 500 m 分 辨 率 NPP 数 据 集 。 

2.2 研究 方法 

2.2.1 技术 路 线 本 研究 以 喀什 噶 尔 绿洲 耕地 为 研 
究 对 象 ,运用 RS/GIS 技术 ,通过 处 理 生长 季 6 期 
landsat 8 OLI 遥 感 数据 提取 植被 指数 3 个 .盐分 指数 
4 人 个、 组 帽 (K-T) 变 换 因 子 2 个 .下 垫 面 反射 因子 2 
个 .组 合 指 数 特征 空间 3 个 波段 反射 率 因 子 6 个 ， 
用 线性 拟 合 的 方法 将 植被 NPP 数 据 进行 降 尺 度 , 并 
通过 土地 利用 数据 提取 研究 区 耕地 开垦 年 限 , 将 20 
个 遥感 参量 同 土壤 采样 及 实测 数据 进行 相关 分 析 ， 
得 到 优选 的 遥感 特征 变量 ;通过 粒子 群 优 化 算法 
(PSO ) 优 化 的 概率 神经 网 络 (PNN ) 模 型 进行 盐 渍 化 
程度 分 类 ,得 到 研究 区 耕地 盐 渍 化 程度 分 布 ,并 将 
耕地 盐 渍 化 程度 与 开垦 年 限 和 植被 生产 力 进 行 琶 
加 ,分 析 耕 地 盐 渍 化 与 耕地 开垦 年 限 的 关系 以 及 耕 
地 盐 涡 化 对 耕地 生产 力 的 影响 程度 ,为 干旱 与 半 干 
旱 区 盐 渍 化 反 演 研究 以 及 当地 盐 渍 化 耕地 管理 提 
供 人 参考 (图 2)。 

2.2.2 遥感 指数 特征 变量 ”候选 的 遥感 指数 特征 变 
量 主要 有 植被 指数 3 个 .盐分 指数 4 个 K-T 变换 因 
子 2 个 、 下 热 面 反射 因子 2 个 组合 指数 特征 空间 3 
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HAE: 基于 PSO-PNN 模型 的 喀什 噶 尔 绿洲 耕地 盐 渍 化 分 析 


土壤 采样 及 实 
测 数 据 土地 利用 数据 NPP(500 m) 


[EUN ORBIS pH 等 


NPP(30 m) 


特征 变量 
选取 
TERE 
征 变量 
盐 渍 化 估 通过 、 
EREA 进行 盐 渍 化 优化 训练 
测 模型 OK d pun 
是 
耕地 盐 渍 化 ; 
VME STR 
程度 分 布 年 限 的 关系 
结果 分 析 | NPP < 又 加 分 析 >| 


耕地 开垦 年 限 


盐 渍 化 对 生 
力 的 影响 力 


TE: NPP 为 植被 净 初 级 生产 力 ;K-T 变换 为 疆 帽 变换 ;PNN 为 概率 神经 网 络 ;PSO 为 粒子 群 优化 算法 。 


图 2 技术 路 线 


Fig. 2 Technical route 


个 波段 反射 率 因 子 6 个 ( 表 1)25 5 : 
223 相关 分 析 为 探究 过 感 特征 变量 与 实地 采样 
点 盐分 之 间 的 相关 性 ,采用 Pearson 相关 系数 表示 各 
变量 之 间 的 相关 性 。 


Ys, -39;-5) 
p z= ml (1) 


| Yes E 


式 中 : r。 为 实地 采样 点 样本 总 盐 含 量 与 遥感 指数 之 
间 的 相关 系数 ; x, 为 实地 采样 点 样本 总 盐 含 量 (mg， 
g); y, 为 实地 采样 点 遥感 指数 ; x 为 实地 采样 点 样 
本 总 盐 含量 平均 值 (mg'.g); 7 为 实地 采样 点 遥感 
BCE LE ; n 为 样本 数 。 
Pearson 相关 系数 的 检验 统计 量 为 : 

_ 

ES 
式 中 :7 为 Pearson 相关 系数 的 检验 统计 量 ;r 为 相关 
系数 ;n 为 样本 数 。 
2.2.4 NPP 数 据 降 尺 度 ”本 文通 过 线性 拟 合 的 方法 
转换 成 30 m 分辩 率 的 NPP 数 据 , 然 后 通过 裁剪 得 到 


T (2) 


研究 区 耕地 NPP 数 据 。 具 体 方法 为 : 

基于 30 m 分 辨 率 的 土地 利用 图 ,将 土地 利用 类 
型 图 根据 需要 重新 分 类 。 首 先 遍历 所 有 像 元 ,判断 
NPP 栅 格 图 像 中 的 每 个 像 元 是 否 是 纯 像 元 ,提取 出 
纯 像 元 NPP 及 其 对 应 的 中 等 分 辩 率 增强 型 植被 指 
数 (EVD) 像 元 值 ;人 然后 将 每 种 土地 利用 类 型 对 应 的 
纯 像 元 EVI 数 据 同 NPP 数 据 进行 线性 拟 合 , 线 性 拟 
合 使 用 的 是 接口 描述 语言 (IDL) 中 的 linfit 函数 ,分 
别 得 到 每 种 土地 利用 类 型 对 应 的 NPP 与 EVI 之 间 的 
WEAR“: 


NPP, =aXEVI, +b (3) 
AP: NPP, Z2 OE SER 28 TH] ABESSE 
777]; EVI, 为 纯 像 元 对 应 的 EVI 值 ; a. b 为 待 拟 合 
参数 。 基 于 拟 合 公式 所 得 的 参数 a、b 及 植被 生长 
季 高 空间 分 辨 率 的 EVI ,实现 研究 区 高 空间 分 状 率 
NPP 的 估算 “*: 


NPP, =aXEVI, +b (4) 
式 中 : NPP, 为 计算 所 得 的 高 分 辩 率 植被 生产 力 数 
据 ; EVI, 为 植被 生长 季 遥 感 影像 经 过 波段 运算 后 得 
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del 遥感 指数 特征 变量 


Tab.1 Characteristic variable of remote sensing index 


类 别 名 称 公式 

植被 指数 归 一 化 植被 指数 (NDVI) (NIR - R)/(NIR + R) 

增强 型 植被 指数 (EVI) has (NIR - R) 

(NIR-6XR-7.5XB+ 1) 

ia a TEM FAC le xNIR+1)- Je x NIR + 1) - 8X (NIR - R)| +2 
盐分 指数 盐分 指数 (SIT) JBXR 

归 一 化 盐分 指数 (NDSI) R-NIR 

R+NIR 
盐分 指数 2(SL) JCE 


K-T 变换 因子 


下 垫 面 反射 


多 维特 征 空间 


2H 


波段 反射 率 


子 


盐分 指数 3(SD) 
湿度 指数 (WI) 
绿 度 指数 (GV1) 

地 表 反 照 率 (Albedo) 
Fe:0; Ett ACC LL) 


MSAVI-WI 特 征 空间 (MWI) 


DVI-SI 特 征 空间 (SDI) 


ULEE2—7(B2,B3.,B4,B5,B6,B7) 


MSAVI-WI-SI 43 (iE 25 [HR] MWSI) 


IG +R +NIR? 

0.1511 x B +0.1972 x G 0.3283 x R + 0.3047 X NIR - 0.7117 x SWIR, - 0.4559 x SWIR, 
—0.16 x B - 0.28 x G- 0.49 x R + 0.79 x NIR + 0.0002 X SWIR, - 0.14 x SWIR, 

0.356 x B + 0.130 x R + 0.373 x NIR + 0.085 x SWIR, + 0.072 x SWIR, — 0.0018 

R 


B 
A(MSAVI - 1f +(WI- 1)? 


(NDVI- I? +SP 


J(MSAVI- 1? +(WI- D +SP 
B. G, R, NIR, SWIR,, SWIR, 


iE K-T AE y ae ; NIR 为 近 红 外 波段 ;R 为 红 波 段 ;B 为 蓝 波段 ;6 为 绿 波段 ; SWIR 为 短波 红外 1 波段 (1.57~1.65 um); SWIR, 为 短波 红 


外 2 波段 (2.11~2.29 pm)。 


到 的 生长 季 30m 分 辨 率 EVI 数 据 ; a. b 为 拟 合 参数 。 
NPP 验证 数据 来 自 于 2000—2019 年 新 疆 统计 
年 鉴 。 其 中 农业 统计 包含 各 县 级 行政 单位 的 数据 ， 
并 具有 连续 记录 ,因此 将 统计 数据 应 用 于 NPP 模 型 
模拟 结果 验证 是 可 行 的 。NPP 验证 方法 是 基于 作物 
产量 和 面积 数据 进行 验证 的 。 该 方法 的 原理 如 下 : 
根据 不 同 作物 收获 部 位 的 含水 量 和 收获 指数 值 ,将 
农业 统计 得 到 的 各 种 作物 产量 转化 为 植被 碳 储量 ， 
AX RU: 
X Y X(1-MCj)x045 /< 
NPP= > es Da (5) 
式 中 :WW 为 作物 类 型 的 数量 ;了 ,为 作物 i 的 总 产量 
(g); 4, 为 作物 i 的 收获 面积 (mx); MC, 为 收获 部 分 作 
物 ; 的 水 分 含量 (gg 9); HI, 为 作物 i 的 收获 指数 。 
经 过 Pearson 相关 性 检验 表明 , 降 尺 度 得 到 的 
NPP 数 据 集 和 了 验证 数据 之 间 的 相关 性 是 显著 的 (R= 
0.64, P«0.001) ,所 以 降 尺 度 得 到 的 NPP 数 据 可 用 于 


本 次 研究 。 
2.3 盐 渍 化 估 测 模型 
2.3.1 概率 神经 网 络 (PNN) 模 型 与 传统 BP 网 络 


相 比 ,PNN 网 络 学 习 简 单 训练 速度 快 、 收 敛 性 好 。 


该 网 络 在 网 络 结构 上 也 更 加 灵活 ,具有 非常 好 的 容 
错 性 ,特别 适合 解决 模式 分 类 问题 。PNN 有 输入 
层 、 隐 含 层 、 求 和 层 和 输出 层 。 通 过 高 斯 函数 运算 ， 
PNN 模 型 的 输出 可 以 接近 非 线 性 。PNN 的 分 类 过 
程 为 :首先 ,需要 将 样本 代入 模型 输入 层 , 然 后 对 样 
本 进行 归 一 化 处 理 以 消除 量 纲 影响 ,并 将 处 理 后 的 
变量 代入 隐 含 层 进行 高 斯 运算 ,进而 将 隐 含 层 的 结 
果 代 入 求 和 层 并 进行 求 和 与 相关 计算 ,最 后 得 到 盐 
渍 化 程度 分 类 并 输出 。 

2.3.2 粒子 群 优化 (PSO) 算 法 平滑 参数 是 PNN 模 
型 训练 的 关键 参数 ,影响 分 类 和 拟 合 预 测 结 果 。 本 
SCRA PSO BIR! 寻找 最 优 参数 ,PSO 算法 的 流 
程 图 如 图 3 所 示 。 

2.3.3 误差 检验 方法 “耕地 盐 渍 化 估 测 模型 一 般 采 
用 交叉 验证 的 方法 与 独立 数据 集 验 证 的 方法 对 佑 
测 结果 进行 严格 的 精度 验证 呈 。 独 立 数据 集 验证 
方法 是 一 种 更 加 直接 的 估计 空间 不 确定 性 的 方法 ， 
验证 的 流程 为 在 所 有 的 样本 中 ,选择 大 部 分 样本 用 
来 进行 机 器 学 习 , 之 后 用 剩余 的 样本 进行 精度 验证 ， 
通过 对 比 模拟 值 与 真实 值 之 间 的 接近 程度 来 判断 估 
测 结果 的 优 劣 性 ,本 文选 择 独立 数据 集 验 证 方法 。 
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计算 每 个 粒子 
的 适应 度 


更 新 个 体 和 全 
局 最 佳 位 置 


图 3 PSO 流 程 图 
Fig.3 Flow chart of particle swarm optimization (PSO) 


3 结果 与 分 析 


3.1 盐 渍 化 估 测 分 析 模 型 建立 

3.1.1 :£ £A € uto 将 采样 点 处 的 遥感 指数 同 总 
盐 进 行 Pearson 相关 分 析 ( 图 4) ,从 每 类 指标 中 选取 
显著 相关 且 代 表 性 强 的 遥感 参量 。 从 网 4 可 以 看 
出 ,在 研究 区 中 ,与 总 盐 呈 正 相 关 的 有 波段 2(B2)、 
波段 3(B3) 波段 4(B4) 波段 6(B6) .波段 7(B7)、 
MSAVI-WI-SI 特 征 空 间 (MWSI) .NDVI-SI 特 征 空 间 
(SDI) , MSAVI- WI 特征 空间 (MWI) 、 地 表 反 射 率 
(Albedo) 、 盐 分 指数 3(SL) .盐分 指数 2(SL) . 归 一 化 
盐分 指数 (NDSI) 、 盐 分 指数 (SD) ; 与 总 盐 呈 负 相 关 
的 有 波段 5(B5) , FejO, 质量 分 数 (J,。) \ 绿 度 指数 
(GVI) 湿度 指数 (WI) 修改 型 土壤 调节 植被 指数 
(MSAVI) 增强 型 植被 指数 (EVI) . 归 一 化 植被 指数 
(NDVI)。 其 中 ,与 总 盐 不 显著 相关 的 指标 有 BS. 


I. o, .MSAVI; f£ P<0.05 水平 下 相关 的 指标 有 SE; 在 
P«0.01 水 平 下 相关 的 指标 有 B3、B2、SDI、 Albedo, 
GVI, SI; 在 P«0.001 水 平 下 相关 的 指标 有 B7、B4、 
MWI、NDSI.NDVI、B6、MWSI、WI、Sb、EVI。 本 文选 
取 植 被 指数 EVI、 盐 分 指数 SL K-T 变换 因子 WI、 多 
维特 征 空间 MWSI 波段 反射 率 因子 B6 共 5 个 遥感 
参量 来 进行 盐 涡 化 程度 反 演 。 
3.1.2 耕地 盐 渍 化 程度 估 测 “在 进行 机 器 学 习 时 ， 
将 EVISL、WIMWSI.B6 共 5 个 遥感 参量 做 为 模型 
的 输入 参量 ,根据 总 盐 将 112 个 样本 平均 分 成 6 级 
( 表 2) , 取 80 个 训练 样本 进行 训练 ,然后 将 剩 下 的 
32 个 样本 用 来 测试 准确 率 。PSO 优化 的 学 习 因 子 
取 1.49445 ,进化 次 数 取 200 ,种群 规模 取 50, 最 终 得 
到 的 最 优 平滑 参数 在 0.5 左右 ,将 平滑 参数 取 0.5 带 
A PNN 预测 模型 中 ,得 到 研究 区 耕地 盐 渍 化 程度 分 
级 , 共 6 级 ,经 测试 集 检验 分 类 正确 率 为 80% 左 右 。 

从 图 $ 可 以 看 出 ,研究 区 内 盐 渍 化 程度 级 别 从 
左 到 右 逐 渐 升 高 。 全 区 盐 渍 化 程度 占 比 最 高 的 是 
第 三 级 和 第 四 级 ,分 别 占 了 研究 区 耕地 面积 的 
37.58% Fil 35.20% ,其 次 为 第 一 级 和 第 六 级 ,分 别 占 
了 人 研究 区 耕地 面积 的 14.78% 和 11.65%, 占 比 最 少 的 
是 第 二 级 和 第 五 级 ,分别 占 了 研究 区 耕地 面积 的 
0.07% 和 0.72%。 除 此 之 外 ,研究 区 耕地 盐 渍 化 水 平 
在 第 一 至 第 三 级 的 区 域 占 比 为 52.43% ,在 第 四 至 第 
六 级 的 区 域 占 比 为 47.37% , 略 小 于 盐 渍 化 水 平 在 第 
一 至 第 三 级 的 区 域 面 积 。 

从 表 3 中 可 以 看 出 ,平均 盐 渍 化 程度 由 低 到 高 
的 县 市 依次 为 阿 克 陶 县 itc ep s FEE ab EL 
英吉 沙县 . 阿 网 什 市 . 岳 普 湖 县 、 伽 师 县 ,其 耕地 盐 
渍 化 水 平 在 第 一 到 第 三 级 的 区 域 面 积 依次 递减 , 耕 
地 盐 渍 化 水 平 在 第 四 到 第 六 级 的 区 域 面 积 依次 递 
增 。 伽 师 县 耕地 面积 最 多 ,有 2162.34 km’, (A FEE 
地 盐 渍 化 水 平 却 不 容 乐观 , 盐 渍 化 水 平 在 第 四 至 第 
六 级 的 区 域 面 积 占 了 60.52% 。 除 此 之 外 , 阿 图 什 
市 、 英 吉 沙 县 .喀什 市 的 耕地 盐 渍 化 等 级 属于 第 六 
级 的 面积 占 比 都 超过 15% ,证 明 在 这 些 区 域 , 尤 其 
是 阿 图 什 市 东部 .英吉 沙县 东部 和 北部 .喀什 市 东 
部 , 盐 渍 化 程度 较 重 。 
3.2 耕地 盐 涡 化 程度 与 开 晨 年 限 和 植被 生产 力 的 
关系 
3.2.1 Se E IG E 57F E 4E X € WHA 
中 可 以 看 出 , 近 10 a 研 究 区 新 开垦 耕地 达 2015.73 
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注 :B7 为 波段 7;B6 为 波段 6;B5 为 波段 5;B4 为 波段 4;B3 为 波段 3;B2 为 波段 2;MSAVI 为 修改 型 土壤 调节 植被 指数 ;SDI 为 NDVI-SI 特 征 空间 ; 
MWI 为 MSAVI-WI 特 征 空 间 ; Io, 为 Fe,0; 质 量 分 数 ;Albedo 为 地 表 反 射 率 ;GVI 为 绿 度 指数 ;WI 为 湿度 指数 ;SE 为 盐分 指数 3;SL 为 盐分 指数 
2;NDSI 为 归 一 化 盐分 指数 ;SI 为 盐分 指数 ;MWSI 为 MSAVI-WI-SI 特 征 空 间 ;EVI 为 增强 型 植被 指数 ;NDVI 为 归 一 化 植被 指数 。 

图 4 候选 指标 与 总 盐 的 相关 性 


Fig.4 Correlation between candidate indices and total salt 


R2 样本 分 级 
Tab.2 Sample grading 

盐 渍 化 程度 等 级 dmg: g” 
第 一 级 0.00-3.05 
第 二 级 3.05-7.42 
第 三 级 7.42~11.30 
第 四 级 11.30~12.60 
第 五 级 12.60-13.90 
第 六 级 13.90~17.60 


km2,45a 以 上 的 老 耕 地 有 4266.72 kms FEX 10 a 新 
开明 耕地 中 ,有 27.03% 的 耕地 盐 渍 化 程度 等 级 在 第 
六 级 , 仅 有 1.27% 的 新 开垦 耕地 在 第 一 级 ,可 见 新 开 
明 耕 地 盐 渍 化 程度 较 高 。 在 耕地 开 星 年 限 大 于 45 a 
的 区 域 ,耕地 盐 渍 化 程度 属于 第 一 级 的 区 域 占 比 为 
22.83% , 仅 有 4.32% 的 耕地 盐 渍 化 程度 属于 第 六 级 ， ”在 东部 ,与 盐 活 化 程度 等 级 分 布 大 概 一 致 de soi 
可 见 耕 地 开垦 年 限 越 大 的 区 域 盐 涡 化 程度 越 低 。 化 程度 等 级 高 的 区 域 ,耕地 开垦 年 限 较 低 。 新 开明 

从 图 6 中 可 以 看 出 ,研究 区 新 增 耕 地 大 都 分 布 ”的 耕地 主要 分 布 在 伽 师 县 . 阿 图 什 市 东部 .后 普 湖 


图 5 喀什 噶 尔 绿洲 耕地 盐 渍 化 分 级 


Fig. 5 Classification of cultivated land salinization in 


Kashgar Oasis 
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表 3 喀什 噶 尔 绿洲 各 市 县 耕地 盐 渍 化 等 级 面积 占 比 


Tab.3 Proportion of the graded area of cultivated land salinization in the cities and counties of Kashgar Oasis 


盐 涡 化 等 级 占 比 /% 
市 县 总 面积 /km 
第 一 级 第 二 级 第 三 级 ”第 一 至 第 三 级 ”第 四 级 第 五 级 第 六 级 ” 第 四 至 第 六 级 
阿 克 陶 县 25.01 0.14 46.38 71.53 21.41 1.67 5.39 28.47 463.77 
Gi 25.05 0.05 42.68 67.78 20.20 1.27 10.76 32.22 1030.86 
喀什 市 22.25 0.00 43.28 65.53 18.07 0.75 15.65 34.47 215.64 
wi 15.40 0.07 43.74 59.21 33.44 0.67 6.69 40.79 1292.31 
英吉 沙县 15.58 0.05 37.43 53.05 27.81 0.86 18.28 46.95 752.85 
阿 图 什 市 20.41 0.00 30.65 51.06 26.89 0.66 21.39 48.94 423.36 
后 普 湖 县 10.92 0.11 33.13 44.16 43.42 0.42 12.00 55.84 985.23 
伽 师 县 6.95 0.07 32.46 39.48 48.53 0.40 11.60 60.52 2162.34 
总 计 14.78 0.07 37.58 52.43 35.20 0.72 11.65 47.57 7326.36 
表 4 喀什 噶 尔 绿洲 耕地 开垦 年 限 和 盐 渍 化 的 关系 
Tab.4 Relationship between the age of cultivated land reclamation and salinization in Kashgar Oasis 
盐 渍 化 等 级 占 比 /% 

SPEED 第 级 第 二 级 第 三 级 第 四 级 第 五 级 第 六 级 
0~10 1.27 0.08 13.24 57.28 1.09 27.03 2015.73 
11~20 3.98 0.04 25.50 53.93 1.09 15.45 486.27 
21-30 9.58 0.11 37.60 42.82 1.29 8.61 167.31 
31-45 12.17 0.05 43.83 34.29 0.83 8.83 390.33 
>45 22.83 0.07 49.89 22.42 0.47 4.32 4266.72 
AMT 14.78 0.07 37.58 35.20 0.72 11.65 7326.36 


图 例 

加 行政 区 界线 
耕地 开垦 年 限 /a 
E >45 

— 30-45 
120-30 
110-20 

mu 0-10 


图 6 MRAP MS AR aD ETT Ek EPR 
Fig. 6 Age of cultivated land in Kashgar Oasis 


FL EB ELAR TB .英吉 沙县 东部 和 北部 。 而 焉 附 县 
西部 、 阿 克 陶 县 新 开展 耕 地 较 少 ,大 都 为 开展 年 限 
在 100 a 以 上 的 老 绿洲 农业 区 。 耕 地 开 居 年限 较 小 
的 盐 溃 化 耕地 主要 位 于 人 研究 区 东部 。 


不 同 地 区 的 气候 .生态 环境 不 同 , 在 新 疆 , 原 本 
就 存在 很 多 盐分 含量 较 高 的 盐 溃 化 土地 ,这 些 地 区 
如 果 不 充 分 进行 盐 溃 化 治理 就 开 呈 ,就 会 自然 而 然 
变 成 盐 汪 化 耕地 。 一 些 自然 环境 和 土壤 质量 较 好 
的 土地 ,如 果 开 居 后 灌溉 用 水 不 合理 也 会 导致 盐 溃 
化 。 人 研究 区 内 的 耕地 开明 年 限 越 小 , 盐 涡 化 程度 越 
重 , 其 自然 因素 主要 是 后 期 开垦 的 耕地 多 位 于 河流 
下 游 和 绿洲 边缘 ,开明 的 多 为 戈壁 滩 之 类 的 未 利用 
地 ,土地 条 件 比较 差 。 早 期 开垦 的 耕地 土壤 条 件 较 
好 ,多 位 于 河流 上 中 游 和 绿洲 内 部 , 开 晨 的 多 为 草 
地 、 林 地 。 人 为 因素 主要 是 随 着 耕地 开垦 年 限 的 增 
加 ,人 们 治理 盐 渍 化 耕地 的 科学 技术 手段 不 断 改 善 ， 
农业 活动 对 耕地 的 改良 和 治理 效果 逐渐 显现 。 另 
外 ,由 于 管理 水 平 不 一 和 自然 条 件 不 同 ,人 研究 区 东西 
部 农业 活动 存在 差异 ,地 势 也 西高 东 低 ,盐分 容易 随 
地 下 水 和 地 表 水 从 西部 流向 东部 ,到 了 东部 平原 地 
区 ,水 流 变 缓 ,土壤 盐分 逐渐 累积 ,使 东部 地 区 更 易 
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发 生 盐 渍 化 ,因此 , 盐 渍 化 程度 较 强 的 区 域 大 都 发 
生 在 研究 区 东部 。 
3.2.2 耕地 盐 汪 化 对 生产 力 的 影响 ”从 表 5 中 可 以 
看 出 ,研究 区 盐 渍 化 程度 等 级 越 低 ,耕地 NPP 越 大 。 
盐 渍 化 程度 等 级 为 第 一 级 和 第 二 级 的 耕地 平均 NPP 
最 大 ,为 718.93 gem ao 盐 渍 化 程度 等 级 为 第 五 级 
和 第 六 级 的 耕地 平均 NPP 较 小 ,为 639.08 gema" 
盐 渍 化 程度 等 级 为 第 三 级 和 第 四 级 的 耕地 面积 最 
大 ,为 5210.91 km ,平均 NPP 为 674.90 g.m?.a。 
从 图 7 中 可 以 看 出 ,研究 区 耕地 NPP 较 高 的 区 
域 大 都 分 布 在 西部 ,耕地 NPP 较 底 的 区 域 大 都 分 布 
在 东部 ,与 盐 渍 化 程度 等 级 分 布 大 致 相反 。NPP 较 
低 的 耕地 主要 分 布 在 伽 师 县 . 阿 图 什 市 东部 .后 普 
WE MERE .英吉 沙县 东部 和 北部 ,而 玻 附 县 
西部 、 阿 克 陶 县 耕地 NPP 较 高 ,因此 耕地 盐 渍 化 对 


表 5 喀什 噶 尔 绿洲 耕地 不 同 盐 渍 化 等 级 
面积 及 对 应 NPP 
Tab.S Area of different salinization grades of cultivated 
land and its corresponding NPP in Kashgar Oasis 


盐 渍 化 等 级 面积 km NPP/g:m® a” 
第 一 级 1063.80 719.33 
第 二 级 4.86 630.39 
第 一 至 第 二 级 1068.66 718.93 
第 三 级 2695.77 696.29 
第 四 级 2515.14 651.98 
第 三 至 第 四 级 5210.91 674.90 
第 五 级 50.13 636.22 
第 六 级 808.47 639.25 
第 五 至 第 六 级 858.60 639.08 


图 7 喀什 噶 尔 绿洲 耕地 植被 NPP 
Fig. 7 Vegetation NPP of cultivated land in Kashgar Oasis 


耕地 农作物 生产 力 的 影响 比较 显著 ,在 盐 渍 化 程度 
等 级 高 的 区 域 ,耕地 NPP 较 低 ,在 盐 渍 化 程度 等 级 
低 的 区 域 , 耕 地 NPP 较 高 。 

$b Oi AEA A BAB BEF, dE REIP BE T TE 
受到 抑制 ,微生物 的 代谢 也 受到 影响 ,从 而 使 土壤 
肥力 变 差 .养分 利用 率 变 低 .土壤 板结 .透气 性 差 ， 
植物 生长 也 受到 抑制 。 另 外 ,盐碱地 中 的 土壤 孔隙 
较 小 ,无 论 是 化 学 性 状 还 是 物理 性 状 都 较 差 ,非常 
不 利于 农作物 的 生长 。 经 实地 采样 观察 ,喀什 噶 尔 
绿洲 绝 大 部 分 地 区 的 耕地 采用 了 滴灌 方式 ,少数 区 
域 采用 渠 灌 ,漫灌 方式 , 盐 涡 化 耕地 中 多 有 少量 杂 
草 ,主要 在 未 耕种 的 区 域 上 ,作物 长 势 较 好 或 盐 渍 
化 严重 的 耕地 中 杂 草 极 少 ,西部 老农 业 区 大 部 分 盐 
渍 化 耕地 中 盐 斑 较 少 或 不 明显 ,植被 覆盖 率 高 。 东 
部 盐 渍 化 耕地 中 盐 斑 明显 ,作物 出 苗 及 长 势 较 差 ， 
往往 还 有 明显 连 片 的 无 植被 覆盖 区 域 。 从 盐 渍 化 
耕地 的 分 布 位置 来 看 ,多 分 布 于 农 区 边缘 .老农 区 
外 部 新 明 区 域 。 另 外 ,有 不 少 地 块 被 弃 耕 ,未 耕种 
的 原因 存在 多 样 化 。 总 体 来 讲 主要 是 当地 年 轻 人 
有 很 多 不 愿意 在 农田 耕作 ,农田 管理 水 平 不 一 ,再 
加 上 当地 气候 干旱 、 降 水 稀少 ,有 伴随 着 风沙 侵袭 
等 自然 灾害 ,很 多 耕地 无 法 继续 耕种 。 总 的 来 说 ， 
耕地 盐 渍 化 是 人 为 作用 下 的 自然 灾害 , 它 对 土地 生 


产 力 有 着 难以 逆转 的 危害 。 
4 结论 


(1) 本 文选 取 EVI、S、WI、MWSI、B6 共 5 个 遥 
感 参量 通过 PSO-PNN 模型 进行 盐 涡 化 程度 反 演 准 
确 率 约 为 80%。 

(2) 耕地 开垦 年 限 越 大 的 区 域 盐 渍 化 程度 越 
低 。 在 盐 渍 化 程度 等 级 高 的 区 域 ,和 耕地 开 虹 年限 较 
低 。 新 开垦 的 耕地 主要 分 布 在 研究 区 东部 , 盐 渍 化 
程度 较 高 。 而 研究 区 西部 新 开明 耕地 较 少 ,大 都 为 
JF EERE 45 a 以 上 的 老 绿洲 农业 区 , 盐 渍 化 程度 
较 低 。 

(3) 耕地 盐 涡 化 对 耕地 农作物 生产 力 的 影响 比 
较 显 著 。 研 究 区 耕地 NPP 较 高 的 区 域 大 都 分 布 在 
西部 ,与 盐 渍 化 程度 等 级 分 布 大 致 相反 ,在 盐 渍 化 
程度 等 级 高 的 区 域 , 耕 地 NPP 较 低 , 盐 渍 化 程度 等 
级 越 低 ,耕地 NPP 越 大 。 

本 文选 取 EVI.SL、WIMWSI.B6 共 5 个 遥感 参 
量 通过 PSO-PNN 模型 进行 盐 渍 化 程度 反 演 ,适用 于 
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研究 区 的 盐 渍 化 反 演 ,但 此 方法 对 于 其 他 区 域 的 盐 
渍 化 反 演 研究 的 普 适 性 还 有 竺 进一步 研究 。 另 外 ， 
本 文选 取 了 Landsat 8 OLI 数 据 的 质量 .影像 的 分 辩 
率 等 会 对 人 研究 结果 造成 一 定 的 影响 , 受 选取 样 点 数 
HE 采样 路 线 通达 性 .前 期 预 判 等 条 件 限 制 ,选取 的 
采样 点 有 一 定 的 局 限 性 ,遥感 影像 处 理 过 程 中 由 于 
影像 质量 等 原因 有 一 定 的 不 确定 性 。 在 以 后 的 盐 
渍 化 反 演 中 可 尝试 不 同 遥 感 影 像 和 不 同 模型 参数 
来 进行 研究 。 
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Abstract: Salinization is one of the main causes of low yields in oasis agriculture and a major constraint and 
barrier to agricultural and sustainable development. In order to improve the productivity of saline-cultivated land 
and to promote the sustainable development of oasis agriculture, and taking the cultivated land of the Kashgar 
Oasis in Xinjiang, China, as the research object, this study used Landsat 8 OLI remote sensing image data to 
extract 20 remote sensing indices. It also calculated the reclamation age of cultivated land in the study area based 
on land use data and downscaled the vegetation net primary productivity (NPP) data simulated by the vegetation 
photosynthesis model using a linear- fitting method. In this analysis, soil sampling and measured data were used 
to obtain the relevant remote sensing characteristic variables, the probabilistic neural network (PNN) model of 
particle swarm optimization (PSO) optimization was used to classify the degree of salinization, and, finally, the 
distribution of the salinization level of cultivated land in the study area was obtained and then superimposed onto 
the cultivated land reclamation age and NPP in the study area. The following conclusions were reached: (1) In 
this paper, five remote sensing parameters: enhanced vegetation index (EVI), salinity index 2 (SL), humidity 
index (WI), MSAVI-WI-SI characteristic space (MWSI), and band 6 (B6, 2.11-2.29 um), were selected to invert 
the degree of salinization using the PSO-PNN model, and this method was found to be effective for salinization 
inversion. (2) The greater the reclamation years of cultivated land, the lower the degree of salinization in the area. 
The newly reclaimed cultivated land is primarily located in the eastern part of the study area. The newly 
reclaimed cultivated land in the western part of the study area is more sparse and is mostly comprised of oasis 
agricultural areas with a land age of more than 45 years. (3) Salinization of cultivated land has significantly 
reduced its productivity. Most of the areas with a higher NPP of cultivated land are located in the west, and most 
of the lower areas are in the east, which is nearly the inverse of the hierarchical distribution of salinization 
degrees. 

Key words: remote sensing (RS); geographic information system (GIS); salinization inversion; particle swarm 


optimization; probabilistic neural network 


